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Abstract 

A unified theory is given of dynamically modified decay and decoherence of field- 
driven multilevel multipartite entangled states that are weakly coupled to zero- 
temperature baths or undergo random phase fluctuations. The theory allows for 
arbitrary local differences in their coupling to the environment. Due to such differ- 
ences, the optimal driving-field modulation to ensure maximal fidelity is found to 
substantially differ from conventional "Bang-Bang" or 7r-phase flips of the single- 
qubit evolution. 
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1 Introduction 



A quantum system may decohere, under the influence of its environment, in 
one (or both) of the following fashions: (a) Its population may decay to a 
continuum or a thermal bath, a process that characterizes spontaneous emis- 
sion of photons by excited atoms (0), vibrational and collisional relaxation of 
trapped ions Q aid cold atoms ta optical lattices §, as well as the relax- 
ation of currenTbiased Josephson junctions gg). (b) It may undergo proper 
dephasing, which randomizes the phases but does not affect the population of 
quantum states, as in the case of phase interrupting collisions (jy). 
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Most theoretical and experimental methods aimed at assessing and controlling 
(suppressing) the effects of decoherence of qubits (any two-level system, that 
is the quantum equivalent of a classical bit) have focussed on one of two partic- 
ular situations: (a) single qubits decohering independently; or (b) many qubits 
collectively perturbed by the same environment. Thus, quantum communica- 
tion protocols based on entangled two-photon states have been studied under 
collective depolarization conditions, namely, identical random fluctuations of 
the polarization for both photons (|7|; |8J). Entangled qubits that reside at the 
same site or in equivalent sites of the system, e.g. atoms in optical lattices, 
have likewise been assumed to undergo identical decoherence. 



For independently decohering qubits, the most powerful approach suggested 
thus far for the suppression of decoherence appears to be the "dynamical 
decoupling" (DD) of the system from the bath fi E3 H E B H H H 



decoupling' 

BEE 




The standard 



'bang-bang" DD, i.e. 7r-phase flips of 
the coupling via strong and sufficiently frequent resonant pulses driving the 
qubit (flit 13t 14), has been proposed for the suppression of proper dephasing 
(|23l). Several extensions have been suggested to further optimize DD under 
proper dephasing, such as multipulse control (flih. continuous DD (18) and 
concatenated DD 0). DD has also been adapted to suppress other types of 
decoherence couplings such as internal state coupling (21) and heating (jl3 ). 



Our grou 



26; 27 



27 



has proposed a universal strategy of approximate DD (J2J; |2f 
29L lict lilh for both decay and proper dephasing, by either 



pulsed or continuous wave (CW) modulation of the system-bath coupling. 
This strategy allows us to tailor the strength and rate of the modulating pulses 
to the spectrum of the bath (or continuum) by means of a simple universal 
formula. In many cases, the standard 7r-phases "bang-bang" is then found to 
be inadequate or non-optimal. 



In the collective decoherence situation, it is possible to single out decoherence- 
free snbspaces (DFS) Q, wherein symmetrically degenerate many-qubit states, 
also known as "dark" or "trapping" states (0), are decoupled from the bath 



17; 33; 34; 3 



Entangled states of two or more particles, wherein each particle travels along 
a different channel or is stored at a different site in the system, may present 
more challenging problems insofar as combatting and controlling decoherence 
effects are concerned: if their channels or sites are differently coupled to the 
environment, is their entanglement more fragile? Is it harder to protect? To 
answer these fundamental questions, we develop a very general treatment. 
The present treatment extends our previously published single-qubit universal 



strategy (j24t l26t 1271 : l36t l37l ) to multiple entangled multilevel systems (parti- 



cles ) which are either coupled to partly correlated (or uncorrelated) zero tem- 
perature baths or undergo locally-varying random dephasing. Furthermore, 
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it applies to any difference between the couplings of individual particles to 
the environment. This difference may range from the large-difference limit of 
completely independent couplings, which can be treated by the single-particle 
dynamical control of decoherence via modulation of the system-bath cou- 
pling, to the opposite zero-difference limit of completely identical couplings, 
allo wing for multi-particle collective behavior and decoherence-free variables 



dl6|; lilt |33t 13J; 135|; l38t l39t |40t |4l() . The general treatment presented here is 
valid anywhere between these two limits and allows us to pose and answer the 
key question: under what conditions, if any, is local control by modulation, 
addressing each particle individually, preferable to global control, which does 
not discriminate between the particles? 

We show that in the realistic scenario, where the particles are differently cou- 
pled to the bath, it is advantageous to locally control each particle by indi- 
vidual modulation, even if such modulation is suboptimal for suppressing the 
decoherence for the single particle. This local modulation allows synchroniz- 
ing the phase-relation between the different modulations and eliminates the 
cross-coupling between the different systems. As a result, it allows us to pre- 
serve the multipartite entanglement and reduces the multipartite decoherence 
problem to the single particle decoherence problem. Throughout the paper we 
show the advantages of local modulation, over global modulation (i.e. identical 
modulation for all systems and levels), as regards the preservation of arbitrary 
initial states, preservation of entanglement and the intriguing possibility of en- 
tanglement increase compared to its initial value. 

In section 2 we present the general formalism, in terms of the systems, their 
couplings to the baths and the modulation used. In section 3 we investigate in 
detail the coupling to zero-temperature baths for multiple multilevel systems, 
and focus on two specific examples, namely a single multilevel system and 
a singly-excited collective entangled-state of many two-level-systems (TLS). 
This is followed in section 4 by a description of multiple TLS (qubits) under- 
going proper dephasing, with a specific example of Bell states. A discussion 
of the results is given in section 5. 



2 General Formalism 



Our total system is composed of M systems, each having a ground state and 
Nj excited states, \g) ■ and |n) ., respectively, where j = 1, ...,M. Each of the 
excited states of each system has a different energy, Uj n . The M systems are 
coupled to a bath and are subject to proper dephasing. Since the coupling to 
the bath may differ from one system to another and for every excited level, 
each is modulated by a different Stark shift H5j tn (t) and a driving field Vj tn (t). 
The total Hamiltonian is the sum of the system (S), bath (B) and interaction 
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(I) Hamiltonians: 



H(t) = H^ s \t) + + HW(t), (1) 

M 

H s {t)=Y,H^t)^ 3 I r , (2) 



j'=i 



n=l 



+hJ2V jtn (t)(\n) jj (g\ + H.c), (3) 

n=l 

H B = hJ2u k \k)(k\, (4) 

M _ Nj 

H\t) =hY,Y,^[hn{t)Vk,3,n{\n){g\) j \vac){k\®I j ^ j 

j=l k n=l 

M Nj 

+H.C] WWiM (5) 

j=l n=l 

Here / is the identity operator, tj tn {t) is the time dependent modulation field 
of the n th excited level of system j, /ik,j,n is the coupling coefficient of the n th 
excited level of system j to the first excited state \k) of the single bosonic 
reservoir mode k, and its proper dephasing Sj n (t) is treated semiclassically. 
H.c. are Hermitian conjugates. The system Hamiltonian includes the system 
terms, as well as the modulation (i.e. Stark-shift and driving field) terms. 
The decoherence effects (both coupling to the bath and the proper dephasing) 
compose the interaction Hamiltonian. 

Two decoherence scenarios will be discussed separately: (i) the coupling to 
the zero-temperature bath is dominant and proper dephasing is negligible, 
e.g. entangled atoms coupled to a cold phonon bath; (ii) proper dephasing 
is dominant and one may neglect the coupling to the bath, e.g. entangled 
photons in coupled fluctuating birefringent fibers (Fig. 1). The treatments of 
both scenarios are based on analogous formalisms. The goal is to optimize the 
modulation or driving in order to ensure maximal fidelity as time goes on. 

The most general state in the system discussed here can be represented in the 
basis of Nt = l\f=i Nj states: 

|*> = £W (6) 
1=1 



In quantum information (QI) implementations it is preferable to use the in- 
teraction representation so that the fidelity of an initial state, |^ ; ), defined 
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as 

F,(f) = |<*,|*(t)>| a , 



(7) 



ensures that in the free (unperturbed) system it remains unity at any time. 




1 2 ... m 



Fig. 1. (a) Entangled multilevel systems with different couplings to a phonon 
bath or different proper dephasings, via <&jj ;nn >(t). Their cross-coupling is through 
Qjj',nn'(t)- The systems are modulated by £j n (t). (b) Several polarization-entangled 
photons propagating through adjacent (coupled) fibers that exhibit fluctuating bire- 
fringence, (c) Entangled systems in tunnel-coupled multi-level wells of a washboard 
potential. There is no direct coupling between the wells at t > 0, only different 
relaxation of each well to the continuum. 

3 Coupling to zero-temperature bath 

3.1 General expressions for dynamical control of zero-temperature decay 

First consider the scenario of different couplings of the systems to a zero- 
temperature bath. The proper dephasing term is neglected in eq. (5), i.e., 
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Sj n (t) = 0, and we set the driving fields to zero, i.e., Vj jTl (t) = 0. 

The difference in the couplings to the k-th mode of the bath is quantified 
by the cross product of their coupling coefficients, (eq. (5)) Hk,j,nHk,j',ri ■ Two 
extreme limits can be discussed: (a) Hk,j,nHk,j' >' — Vk,j ^ j',n,n', in which 
the sets of {kj} modes are separately coupled to each system, making the total 
hamiltonian separable into contributions of the M systems; (b) Hk,j,n = l^k,j',n' 
Vfc, j, f ,n,n' , meaning that the systems are identically coupled to the bath. 

There is initially one excitation in our total system, thus the full wave function 
in this scenario is: 

M M Nj 

Mi)) = E«S(*)I*> ® 1^- + EE a i>B (f)|n»oc) ® \g) f (8) 
k j=l J=l n=l j'+i 



where \vac) is the vacuum state of the bath. In order to analyze the time- 
evolution of the wave function, written as a column vector a(t) = {cxj >n (t)}, 
it is expedient to express it in the interaction picture, 

a jtn (t) = e^-w/o^Wo,,^. (9) 



The Schrodinger equation for the coupled {(Xj,n{t)} and {«($(£)} amplitudes 
in (6) may be reduced, upon eliminating the {«o(^)} amplitudes and trans- 
forming to the interaction picture, to the following exact integro-differential 
equation: 

t 

&iM =J dt 'T, " t')Kf,nn'(t, t'y^^'^'a^At') (fO) 

3',n' 



Here (t) is the reservoir- response matrix, given by: 

^,nn'(t) = J duGfwMe-*", (If) 

Gfj', nn >(u) = ^~ 2 ^2^kd,nfJlj', n 'S(u - u k ). (12) 
k 

and K D (£,£') is the modulation matrix, given by: 
where 

4n(t)=hn(t)e- l f> TS ^ T) (14) 
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accounts for both the modulation and the Stark shift. 

Since we are interested in modulations yielding slowly varying solutions, we 
can pull aj' >n r(t') ~ ajiy(t) out of the integrand in eq. (10), thereby reducing 
it to the differential equation: 

& = -W D (t)dt. (15) 

Here the dynamically-controlled decoherence matrix, W D (t) = {Wjj, nn ,(t)}, 
is a convolution of the modulation and response matrices: 

t 

W°, nn ,(t) = J dt'^ nn ,(t - t')K°, nn/ (t, t')e^-^'^ (16) 
o 

The solution to eq. (15) is given by: 

a(t) = T +e - jD( *' ) a(0) (17) 

where J D = {Jf r , nn ,} and 
t t< 

J?f,nn'(t) =Jdt'J dt"^, nn ,(t' - t")K°, nn ,(t', O^"^ (18) 


Equations (15)-(18) are the most general expressions possible for the deco- 
herence of multilevel, multipartite entangled state under dynamical control by 
modulation. In what follows we explore the consequences of these expressions. 

Two alternative control strategies may be conceived of. The first one is that 
of global modulation, meaning the modulation is identical for all systems, 
e f,n = e f',n' ^j-ij'i n i n ' ■ m this case, the decoherence matrix (16) retains its off- 
diagonal elements and the different states mix. The alternative strategy is that 
of local modulations, i.e. e® n ^ e®, n , Vj, f, n, n'. It will be shown advantageous 
to equalize the rates of decay of all systems {j} and all levels {n}, and to avoid 
their mixing by the decoherence. These requirements amount to fulfilling the 
following conditions: 



JS',uu' = Vj^j'orn^n' (19) 
exp[-Jj? jlin , n ,{t) - % Jq dt'5 fin/ (t')} 



eM-J^nn^-iJodt'SjAt')} 



= 1 Vj,f,n,ri (20) 



33 

which means that 
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»j£nn(*) = 3M5,ii Vj,n (21) 
i < 

%J?j,nn(t) + J dt'5j, n (t') = 3Jg >u (f) + / [mod 27T] (22) 



These conditions imply that different modulations must be applied to each 
system, in all cases, whether the systems are coupled to the same bath or to 
different baths. Our ability to fulfil these conditions and, at the same time, 
minimize the decay /decoherence of amplitudes aj jTl (t), can be quantified in 
terms of the mixing Cj >n and decay A(t) parameters: 

Cj, n (t) = a hn {t)/a hl {t) (23) 
A(t) = a 1:1 (t) lY,\c jn (t)\ 2 - (24) 

If only condition (19) is met, then eq. (23) yields 

In what follows we distinguish between two possible objectives: (i) the preser- 
vation of the initial multipartite entangled state; (ii) the steering of a partly- 
entangled (or unentangled) initial multipartite state to a fully multipartite 
entangled state, both in the presence of decoherence and modulation. 



3.2 Preservation of an initial entangled state 



If one wishes to preserve an initial entangled state, then one should impose 
condition (20), whereby the different states do not mix, i.e. Cj jn (t) = Cj 5 „(0), 
but rather decay at a modified rate, J D (t) = J® tll (t), where 

|tf(t)) = e - jD (*)-*/o*'«i,i(f)|^(o)) (26) 



The fidelity under these conditions is identical for all initial states and is given 
by: 

Fi(t) = e' mjD V (27) 



Expression (26), obtained under conditions (19)-(20), is our result for the 
optimal fidelity of preservation (27) under zero-temperature decay: namely, 
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optimal preservation requires the elimination of state mixing and equal sup- 
pression of the decay for all systems. 



3. 3 Steering 



If one wishes to steer an initial state |\l/(0)), to a desired state \^ d ), it is 
possible to exploit the local modulation and the different decoherence rates in 
order to acquire the desired state at a specific time t. The resultant fidelity 
is then defined as F d (t) = |(\l/ d |\l/(t))| 2 . In order to effectively control the 
amplitude ratios of the states, while avoiding their undesired mixing, it is 
expedient to define the mixing parameters as 

c jt „(t) = a,>(t)/ao» max (t) (28) 



where (j, n) max is chosen such that «(j,n) max is the largest amplitude of the 
desired state. This choice ensures that each system j and each level n are con- 
trolled independently (locally) without affecting the other systems. If (j, n) max = 
(1, 1), then condition (19) yields 

exp[-jg ;11 (t)-^ ^tV(tO] _ 



exp[-Jg, rm (t)-iJ*dtf8 jtn (t')] ~ c d hn ' 



and at time t the state has evolved to: 



|tf(t)) = e -J D W-«/o^i.i(O ai - 1 (°) |^ 



) (30) 



where J D (t) = Ju,u(t) is given in eq. (18). The resulting fidelity is 

F d (t) = e-™ jD V\a hl (0)/<i\ 2 - (31) 



Expression (30), obtained under conditions (19), (29), yields our result for the 
optimal steering towards entanglement fidelity under zero-temperature decay: 
matching the decay rates and initial mixing parameters to the desired mixing 
parameters. 

Next, the implications of the foregoing general recipes for state fidelity preser- 
vation and steering will be analyzed for the following scenarios: 
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3.4 Example: A single multilevel system 



In this case there is a single system with N excited levels, thus the j subscript 
will be omitted. The bath response matrix will be taken to be: 

*°M = Cr m <d n .e-*' 4 *° (32) 



where c nn i is a constant coupling matrix d n = cosr] nj with r] n being the angle 
of transition dipole, and t c is the correlation time of the bath. Here, the system 
eigenstates are equidistant, u n = iO\ + (n — 1) A. We shall use impulsive phase 
modulation, e^(t) = e^ t ' T ' n . Here [...] denote the integer part, r is the pulse 
duration and 9 n is the phase change of level n. 

In figure 2 one can observe the symmetrization of J® n i{t) as a function of 
time. The system has N = 4 levels. By choosing 9 n such that the long-time 
limit of J® n is the same for all levels, one achieves the elimination of the 
off-diagonal terms by different modulations, and the symmetrization of the 
diagonal elements. 




10 20 30 

t 



Fig. 2. J as a function of time. Here lo\ = 0.5, A = 0.1, t c = 1, 
Vn/^ = {0.246,0.0,0.326,0.370}, and c nn = 1.0, c nn i = 0.5. The modulation in- 
terval time is r = t c and the phase changes are 9 n /ir = {1.0, 9.0, 8.0, 7.0}. The blue 
lines denote the 4 diagonal elements of J D , while the red lines are the off-diagonal 
ones. 

Figure 3 displays the decay and mixing parameters as a function of the power 
W n /r invested in the impulse phase modulation. For global modulation, 6 n = 
it is the optimized modulation phase for each level coupled to a Gaussian 
bath, whereas for local modulation, the phase modulation for each level 6 n at 
a given r is found such that symmetry is achieved, when possible. Due to the 
simplicity of the modulation and the large differences in the coupling to the 
bath, symmetrization was not always possible. The x-axis units are those of the 
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mean 9/r = (1/N) J2n@n/i~. As can be seen, one does not increase the decay 
by using local modulation compared to global modulation. However, whenever 
symmetrization is possible, local modulation achieves greater preservation. 

(a) 

1, , , , , , -, 



< 0.5- 



2 4 6 8 10 12 

9/x 

(b) (c) (d) 



1 

^0.8 
0.6 

5 10 5 10 5 10 

dtx e/x e/i 

Fig. 3. Decay parameter \A\ (a) and mixing parameters (02,3,4! (b,c,d respectively) 
as a function of power invested in impulse phase modulation. For global modulation 
(dotted) 9 = 7T whereas for local modulation (solid) for each r the modulation phase 
6 n is chosen such that symmetrization is achieved (whenever possible). Here t = 50 
and other parameters are as in figure 2. 

3.5 Example: Entangled states of M qubits 

Consider a system composed of M two level systems (TLS) or qubits, with 
ground and excited states, \g) and \e), respectively. Only single-excitation 
states are considered here, so the full wave function is given by : 

M M 

m=®\g) ] +Y.di\D? 1 ) (33) 

j=l 1=1 

M 

lA M ) = E^ ) ^|e) J (g)|^ (34) 
i =1 j'¥=j 

where [Df*) provide a (completely entangled) basis for all possible single- 
excitation states, qf = e 2m ^ l ~ l ^ M /yfM. \D^) are zero sum amplitude states, 
except for \D^), which is the symmetric (Dicke) state. We assume here that 
the TLS have the same excitation energies, hu . 

Using the definitions in eqs. (23)-(24) of the decay and mixing parameters, 
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the fidelity of an initial entangled-state is a product of two terms, one due to 
decay and the other due to mixing with other basis states: 

\T M a {l) c (t)\ 2 



One can maximize the fidelity by either reducing the decay, or reducing (or 
even eliminating) the mixing of the entangled-states. For the latter, one has 
to diagonalize and equalize the diagonal elements of the decoherence matrix, 
fulfilling the conditions in eqs. (19)-(20). If these conditions are met, the singly- 
excited entangled-states do not mix with each other. 

In the following numerical example, M = 3 TLS are initially prepared in a 
symmetric singly-excited Dicke state and are coupled to a zero-temperature 

_ t 2 /4t 2 -t 2 /4t 2 , 

bath. The bath response matrix is taken to be &jj>(t) = 7- — — — where 
7 is a coupling constant, tj is the correlation time of TLS i, tq is an arbitrary 
distance and tjj> = \Tj — ry|, where rj is the position of system j. This model 
may describe residual absorption and scattering (out of their initial modes) 
of three polarization-entangled photons in adjacent nearly-overlapping fibers 
(Jiil ) or entangled, vibrationally relaxing atoms at three inequivalent adjacent 
traps or lattice sites, all coupled to the same continuum (fig. 1) (ji^ l. Impulsive 
phase modulation is used as before, ef(t) = e % ^ T ^ 3 . 

Figure 4 shows the mixing and decay parameters (eq. (35)) as a function 
of time, for three TLS with cross-coupling between their relaxations, with 
Yj = {r cos(27rj'/M), r sin(27rj/M), 0.0}. For an identical coupling of all 
qubits (TLS) to the bath, i.e. equal correlation times tj = tjr, one sees that a 
global modulation, meaning the same modulation for the three TLS, results 
in zero mixing, whereas local, or different, modulation results in increasing 
mixing with time. However, for the case of different couplings to the bath, i.e. 
tj 7^ tjr, local modulation can eliminate the mixing, if we choose the optimal 
modulation that equalizes the diagonal decoherence matrix elements, whereas 
global modulation results in an increased mixing with time. 

For any difference in the coupling of the qubits, the results are qualitatively 
similar for all initial singly-excited entangled-states {D^ 1 ). Thus, for identi- 
cal couplings, local modulation achieves similar decay with increased mixing, 
whereas for different couplings, local modulation reduces both decay and mix- 
ing compared to the known global modulation (th e so-called "parity-kicks", 
i.e. 7r-phase flips for all four qubits|To1; Hll : 14; Eh). The optimal recipe is, 



then, to apply M synchronous pulse sequences to the M qubits, but with 
locally- adapted pulse areas: Oj/n = {1.0,0.70,0.58} in the example of fig. 4. 

Figure 5 shows the decay and mixing parameters of steering an initial state 
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20 40 60 80 20 40 60 80 



t t 

Fig. 4. Entangled-state preservation, (a) Decay parameter |^4(t)| as a function of 
time. (b,c) Mixing parameters |c2,3 (i) | as a function of time. The response ma- 
trix parameters are 7 = 0.05, tq = 1.0, the system parameter is u = 0.5 and 
the modulation parameter is r = 1.0. The dotted (solid) lines indicate global (lo- 
cal) modulation, with Oj/tt = 1.0 (fy/vr = {1.0,0.70,0.58}) . The blue (red) lines 
are for identical (different) coupling to the baths for the two TLS, with tj = 1.0 
(i,- = {0.75,0.81,1.0}). 

which is a superposition of several (completely entangled) basis states, with 
c i(0) = {1.0, 1.57, 1.64} and A(0) = 1.0. The desired state is taken to be the 
symmetric singly-excited Dicke state. One can see that using local modulation 
with different pulse rates on the three qubits causes the mixing parameters to 
approach their desired value (i.e. c d 2,3 = 1.0). This result shows that one can 
exploit the different decoherence of the qubits in order to steer a general initial 
state into any desired state. 



4 Proper dephasing 

4-1 Fidelity of M-qubit states 

Consider the scenario where the proper dephasing terms in eqs. (3), (5) are 
dominant, the coupling to the bath and the Stark shifts being neglected. Here 
we assume from the outset to have M TLS (qubits) which undergo different 
proper dephasings. The driving fields of each qubit, Vj(t) are used to dynam- 
ically reduce the proper dephasing. 

In this case, the Hamiltonian is separable into parts pertaining to the individ- 
ual qubits. The wave function of any of them is given by: 

\*)i=W)\9)i + Mt)\e)j (36) 
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Fig. 5. Steering towards a symmetric Dicke state, (a) Decay parameter \A(t)\ as 
a function of time. (b,c) Mixing parameters 1 02,3 (t) | as a function of time. The 
response matrix parameters are 7 = 0.05, tq = 1.0, tj = {0.75, 0.81, 1.0}, the system 
parameter is u = 0.5 and the modulation parameter is r = 1.0. The dotted (solid) 
lines indicate global (local) modulation, with Oj/w = 1.0 {Oj/n = {0.80, 0.56, 0.47}). 



Assuming, for simplicity, that the driving fields are resonant with real enve- 
lope, i.e. Vj(t) = Vj (t)e~ lu ' it + c.c, one can change to the rotating frame 
by defining f3j e {t) = e ZUJ:it Pj e (t) and $j g (t) = /3j 9 (t) and use the rotating wave 
approximation, neglecting terms oscillating at optical frequencies. The diago- 
nalizing basis of the system Hamiltonian of each TLS is: 



t), 



l To {e-^n + 19),) , I I), = ^ (e-^H - l*>i) 



V2 



(37) 



Each TLS wave function is now described by 

l*> i = ^+(*)IT> J - + /3i-(*)U> i . 



(3f 



where 



(39) 



which results in the single-TLS dynamical equation: 

Pi±(t) = TiVj°\t)(3 j± (t) - (& + (f) + &_(t)) 



(40) 



The full wave function is composed of N T = 2 M basis states, I = l...N T , 
which can be presented in a binary representation, meaning I = b\b l 2 ...b l N , with 
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b l j — 0, 1. Here zero denotes f or plus sign and one denotes j or minus sign. 
Thus each basis state |\E^) is a product of M TLS states: 



M 

id' 



W = ®\% (41) 



and the full wave function is given by: 

N T 

l*> = £A|**> ( 42 ) 



where 



M 

a = n p& (43) 



In order to solve for the wave function (6), it is useful to define the column 
vector /3 = {A} and adopt the matrix formulation. We next transform the 
column vector to account for the driving fields: 



P = e iJ dt>P(t>)p (44) 

M 

Pii'(t)=6 w ^2(2b l j -l)V j (t) (45) 

3=1 

where 5u> is Kronecker's delta. This vector fulfills the following dynamical 
equation: 

= -(t/2)W p (t)P (46) 



The transformed proper-dephasing matrix can be split into a part that is 
proportional to the identity I and an off diagonal part: 



(47) 
(48) 
(49) 



Here \l — l'\b is the binary distance between / and I', measuring how many 
qubit-flips are required to get from I' to /, jV->/ is the qubit required to flip 



W p (t) = 


M 

iE^) + w^(t), 

3=1 


<; off (t) = 




0,(t) = 


t 

2 J dt'vj°\t') 
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in order to get from V to / and the sign function sf~*' = b\ — b l - is +1 
for a qubit flip 1 — > and —1 for a qubit flip — > 1. A qub it-flip of qubit j 
means a change of | |) ■ <-> | J,) ■ (eq. (37)) and not the information qubit flip 

l e >i ~ 

The off-diagonal terms of the transformed proper-dephasing matrix are non- 
zero only for elements which require a single qubit flip, and are equal to the 
product of the proper-dephasing rate of the qubit flipped and the modulation 
phase of that qubit, with the appropriate sign e^^'^Sj ,(t). 

Since the proper dephasing term is stochastic, one must define the first and 
second ensemble-averaged-moments, as Sj(t) = and ®fj>(t) = <5j(t)<5j,(0), 

respectively, and adapt the solution to the density matrix p(t) = j3(t)(3^(t). 
Solution of eq. (46) to second order in <5J then corresponds to: 

p(t) = p(0) - - J dt'j dt"[Wr f (t% [Wr ff (t'% P(0)]] (50) 
o o 



It describes the evolution of the density matrix under two consecutive (virtual) 
qubit flips, i.e. excitation and deexcitation, ending up with the same number 
of excitations, but with a random (stochastic) phase. 

The fidelity of an initial basis state, Pfc(0) = l^)^! is now defined as: 

F t (t) = (nV(t)\%) (51) 



It is identical for all initial basis states and is found to be: 



i M 

^(*) = i-«E*^(*) 

t t' 



J?,(t) = Jdt' J dt"^ r {t' -t")Kf r {t\t") 



(52) 

(53) 
(54) 



where = e l ^ j ^\ Since the basis states are product states, this fidelity does 
not pertain to entanglement. Different modulations of individual qubits do not 
affect it. In order to explore the implications of local and global modulations, 
we shall revert to the basis of entangled Bell states. 
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4-2 Fidelity control of Bell states under proper dephasing 



4-2.1 General recipe 

Let us take two TLS, or qubits, which are initially prepared in a Bell state. 
We wish to obtain the conditions that will preserve it. In order to do that, we 
revert to the Bell basis, which is given by 



\B^) = l/V2e^ t (\e) 1 \g) 2 ±\g) 1 \e) 2 ) 
\B, A ) = l/^2(e i ^\e) 1 \e) 2 ±\g) 1 \g) 2 ). 



(55) 
(56) 



This is done by applying the proper rotation matrix to eq. (50). For an initial 
Bell-state p ; (0) = \Bi){Bi\, where I = 1...4, one can then obtain the fidelity, 
F,(t) = (S,|p J (i)|S,),as: 



Ftf) = cos(0±(t))» 

t f 



33 



4>M = J dt 'J dt'% l (t'-t")Kf J , l (t',t") 
o o 

K? jtl (t,t>) = ef*(t)ef(t') 
Kf r ^t') = -Kf r ^t') = e^{t)4\t') 

^(t,o = -^, 2 (t,o = cf(0^(O 



(57) 

(58) 

(59) 
(60) 
(61) 



where 4>±{t) = (<f>i(t) ±<f>2(t))/2 and the 0+ corresponds to k — 1, 3 and 0_ to 
k = 2,4. 

Expressions (57)-(61) provide our recipe for minimizing the Bell-state fidelity 
losses. They hold for any dephasing time-correlations and arbitrary modula- 
tion. 



4-2.2 Numerical example 

In the next numerical example, the response matrix is taken to be 



$£.,(*) = 7e ^ 2 V «' 



(62) 



where 7 is a coupling constant, tj is the correlation time of TLS i, and r^y = 
\tj — Tji\, where Tj is the position of particle j. This model may again describe 
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multi-photon (6) or multi-atom decoherence, as above (fig. 1). Impulsive phase 
modulation is used as before cf(t) = e*'*' 7 ^. 




10 20 30 40 50 
t 




Fig. 6. (a) Fidelity as a function of time. The blue (red) lines indicate triplet (singlet) 
initial states, whereas the solid (dashed) lines show the effects of global (local) 
modulation, with 62 = 0.9tt {62 = O.871"). (b,c) Dephasing rates of singlet (b) and 
triplet (c) initial Bell-states, as a function of modulation phases, B\ % 2- The response 
matrix parameters are 7 = 0.01, t\ = £2 = 1-0 and tq = 1.0 and the modulation 
parameters are t± = T2 = 1.0, B\ = 0.97T. 



Figure 6 shows the fidelity, Fi(t) as a function of time for the singlet, I-B2) and 
triplet, I-B4) states under global and local modulation. The global modulation 
is seen to affect the different states in a different manner: the singlet state 
decoheres more slowly than the triplet state. However, the local (different) 
modulation for the different TLS, eliminates the cross-coupling terms and 
equalizes the decoherence rates of the two states. One can further see, by 
inserting the aforementioned response matrix into (58), that local modulation 
has the same effect as the decorrelation of the two TLS-bath couplings, i.e. 
each of the entangled qubits now decoheres independently, as for r 12 — > 00. 

Thus, locally induced decorrelation of the dephasings can either reduce or 
enhance the Bell-state fidelity compared to standard global ("Bang-Bang") 
7r-phase flips, depending whether the correlated dephasings interfere construc- 
tively (for triplets) or destructively (for singlets). 

If we use singlet and triplet intermittently to encode information, it is advan- 
tageous to use local modulation (in fig. 6: 0i = 0.97T, 9 2 = O.Sn) to equate their 
fidelities, rather than the standard "Bang-Bang" . 
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5 Conclusions 

In this paper we have expounded our comprehensive approach to the dynami- 
cal control of decay and decoherence. Our analysis of multiple field-driven mul- 
tilevel systems which are coupled to partly-correlated or independent baths or 
undergo locally- varying random dephasing has resulted in the universal formu- 
lae (15)-(16) for coupling to zero-temperature bath and (48)-(50) for proper 
dephasing. The merits of local vs. global modulations were presented and are 
summarized below: 

• For different couplings to a zero-temperature bath, one can better preserve 
any initial state by using local modulation which can reduce the decay as 
well as the mixing with other states, more than global modulation. For a sin- 
gle multilevel system, it was shown that local modulation which eliminates 
the cross-decoherence terms, increases the fidelity more than the global mod- 
ulation alternative. For two TLS, it was shown that local modulation better 
preserves an initial Bell-state, whether a singlet or a triplet, compared to 
global 7r-phase "parity kicks" . 

• One can exploit the different couplings to a zero-temperature bath and 
local modulation in order to steer an initial partly-entangled or unentangled 
state, to a desired entangled multipartite state. The ability to match the 
decoherence rates to the desired mixing parameters is made possible by 
using local modulation, which results in lower fidelity losses compared to 
global modulation. 

• Local modulation can effectively decorrelate the different proper dephasings 
of the multiple TLS, resulting in equal dephasing rates for all states. For two 
TLS, we have shown that the singlet and triplet Bell-states acquire the same 
dynamically-modified dephasing rate. This should be beneficial compared 
to the standard global "Bang-Bang" (7r-phase flips) if both states are used 
(intermittently) for information transmission or storage. 

Our general analysis allows one to come up with an optimal choice between 
global and local control, based on the observation that the maximal suppres- 
sion of decoherence is not necessarily the best one. Instead, we demand an 
optimal phase-relation between different, but synchronous local modulations 
of each particle. 
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